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Abstract. 

In this paper the concept of tomography of a coUisionless stellar system of general shape 
is introduced, and a generalization of the Projected Virial Theorem is obtained. Applying the 
tomographic procedure we then derive a new family of virial equations which coincides with the 
already known ones for spherically symmetric systems. This result is obtained without any use 
of explicit expressions for the line-of-sight velocity dispersion, or spherical coordinate system. 
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A very important problem of stellar dynamics is the mass determination of stellar 
systems using their line-of-sight (l.o.s.) velocity dispersion af^^ and surface bright- 
ness profile I. A useful mass estimator, based on the second order Projected Virial 
Theorem (PVT) for spherical coUisionless stellar systems, is obtained integrating 
over the whole projection plane af^^ using I as weight function, (cf. e.g., Ciotti 
and Pellegrini, 1992). Indeed the classical PVT (see, e.g., Kent, 1990) states that 
for any spherical, stationary, non-rotating system with constant mass-to-light ratio 



where R is the projected radius. In (1) the potential energy W = —GM^/rc is 
expressed using the gravitational radius of the system, re, and the last identity is 
just the second order Scalar Virial Theorem, where the kinetic energy is written 
using the (three-dimensional) virial velocity dispersion, cry. In words, the integral 
of the l.o.s. velocity dispersion, weighted on the surface brightness profile of a 
spherical non-rotating galaxy with constant mass-to-light ratio, is equal to 1/3 of 
its (three-dimensional) virial velocity dispersion. The proof of (1) can be obtained 
using for Icr^^g the integral expression given in equation (4.61) of Binney and 
Tremaine (1987) (hereafter BT87, p. 208), changing the order of integration, and 
integrating by parts using the Jeans' equation in spherical coordinates, as shown for 
example by Tonry (1983). The many steps of the proof can obscure the simplicity 
of the result, whose importance is intuitively clear: the potential energy is related 
to the total mass of the galaxy via the knowledge of re, which in principle, for a 
one-component galaxy (i.e. with no dark matter), can be derived by deprojecting 
the surface brightness profile. Because of its importance, it is interesting to extend 
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this result to systems of general shape, in order to obtain a major insight in the 
subject and to generalize (1). For instance Kent (1990) and Merrifield and Kent 
(1990) derived higher order PVT for spherical systems, but the calculations are 
increasingly complicate. 

In this paper a general approach to this problem is presented introducing the 
concept of tomography of a coUisionless stellar systems of unspecified shape. In 
Section 2 the main properties of the projection operator are elucidated. In Section 
3 the iV-th order PVT for systems of any shape is obtained, and in Section 4 
the concept of tomography is presented, together with a theorem on the angular 
mean of the iV-th power of the velocity component along the line-of-sight direction. 
Finally in Section 5 the tomography is applied to the generalized PVT's; in the 
particular case of spherical symmetry the family of any order PVT's [including (1)] 
is recovered, without using explicit expressions for the l.o.s. velocity dispersion, or 
spherical coordinates. 



2 The projection operator 

Let us consider a one-component coUisionless stellar system, described in an orthog- 
onal reference system S = (0;ei,e2,e3) by its density /?(x; f) = £z/(x;f), where 
1/ is the density of light. In S we introduce the orthogonal observer's system 
5" = (O'; fi, £2, fa), and we assume that O = O' at any time, and that S and 5" 
are inertial. Throughout the paper the repeated indices sum convention is used, 
so that X vector in S and ^ = ^jfj is the same vector in 5". Following 

Binney (1985) we do not use the usual Euler angles to specify the orientation of 
5" with respect to S but instead we apply a 3-2-3 rotation^. In this way 1? and (p 
coincide with the polar coordinates of £3 in S, considering to be the polar axis 
and ei the azimuthal one: 

i = Rx, (2) 

with R = R3(^)R2(i?)R3(¥'). Then the unit vector n = £3 from O to the observer 
is given by 

= (cosy? sin 1?, sin sin 1?, cosi?). (3) 

Every observed property Fios associated with the light distribution of a stellar 
system is the projection onto the plane tt' = (^i,^2)-Ln of i^(x; f) (see Sect. 3), 
weighted with the system light density: 

^(^1,6; 0^/o«(^i, 6; t)= [ i/(R^^; t)F{K^e, t) d^3. (4) 

R being orthogonal, the projection operator satisfies the following relation: 

/ I{$i,$2;t)Fios{$u$2;t)d$id$2= [ !/(x; f)^(x; f) rf'x. (5) 

^ Starting from two coincident systems, the Euler angles are defined by a rotation around 83 
{f), then around (■!?) and finally around 83 (■^). In this 3-1-3 rotation the angles ip and are 
not the usual spherical angular coordinates of 83 in S. 
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3 



3 The Projected Virial Theorems 

The definitions introduced up to now are all geometrical, i.e. they do not involve 
the dynamics of the system. All the dynamical properties of a coUisionless stellar 
system are known if its distribution function / = /(x, v; i) is known. This function 
obeys the CoUisionless Boltzmann Equation (see, e.g., BT87, p. 192): 

^ + ^^.^_^^ = 0. (6) 
dt ' dxi dxi dvi 

In fact every macroscopic function -F(x; t) (i.e. describing some property of a 
small element of the system but containing a large number of particles) is defined 
by integrating on the velocity space its microscopic counterpart F[v), using / as 
weight function: 

p{^;t)F{^;t)= [ /(x,v;f)i^(v)dV (7) 
Integrating (6) over 3?'^ after multiplication by F, and considering (7) one obtains: 



dpF^dpv^^ fdF\dl^^_ 
dt dxi \dviJ dxi 

Since we are looking for a generalization of the classical PVT, where the second 
order monomials made with the velocity components are used, we assume for F 
the iV-th power of the v component along n: 

F{^r)=<n,^r>^=v^, (9) 

where < , > is the inner product in the euclidean space 3?'^, and iV > is an 
integer. For example, for N = 1 

t;(n,x;f) = < n, v > =< n, v(x;f) >, (10) 

and for N = 2, defining the second order velocity dispersion tensor to be afj = 

v^n, x; t) = v\n, x; t) + a'p{n, x; t), (11) 
where (t^ = ninjafj. The expansion of (9) for iV > is: 



^\k)\l)^^ ^ "2)(^3 ^1 ^2)- (12) 
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It is evident that the addenda in (and so in v^) are the products between 
the l.o.s. independent velocity monomials and the angular coefficients that on the 
contrary depend only on n. Projecting one obtains 

/(^i,6;iK/o.(^i,6;i)= / ^^^^^6, (13) 

and integrating the l.h.s. of (13) over tt', one finds from (5) 

/ /(^i, 6; t)^ios{^i,^2; t) dir = [ p^d""^ = 2Kn,Ios (14) 

where Kjm^Ios is the projected kinetic energy of order N , which depends on the 
direction of the l.o.s. through the n angular coordinates. For example, for N = 1 
this quantity is the projection of the total momentum (and is according to the 
total momentum conservation), and for iV = 2 it is the usual projected kinetic 
energy. From (8), the equation for which is necessary to derive the generalized 
PVT is: 

^^+^^^^^ + jV„.,?^|*=0. (15) 

Ot OXi OXi 

In the specific case of a stationary system, after multiplication of (15) by Xk and 
integration on 3?'^ one obtains three equations. Multiplying each of these equations 
for rifc, and summing, we finally obtain the iV-th order PVT: 

Theorem For N >0, 

2Kn+i^Ios = -Nn,nkWN-i,zk = -NWn-i,Ios (16) 

where 

WN-l,^k = - [ p'^Xk^d^l^ = Wn-1M- (17) 
JSR3 dxi 

For iV = 1, the r.h.s of (16) do not depends on velocity, Wo^ik being the usual 
potential energy tensor (see, e.g., BT87, p. 67), and one obtains the classical 
second order PVT. 



4 The angular mean operator 

Let us suppose that the observer can move all around the system, observing some 
projected property Fios that generally will depend on the l.o.s. direction. This 
dependence can be eliminated calculating the angular mean of Fios over the solid 
angle: 

\Fiosh = T- / Flos d'n, (18) 
47r 

where now [-F;os]n only depends on the intrinsic (3- dimensional) properties of the 
system. We will call this procedure tomography. For spherically symmetric system 
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the result clearly coincides with the projected property along any l.o.s. direction. 
Then if we succeeded in the explicit calculation of (18), we have at the same time 
(a) a generalization for systems of unspecified shape of the known results for spher- 
ically symmetric systems and (b) a proof of the relations holding for spherically 
symmetric systems avoiding use of particular algebraic expressions between the 
quantities to be projected (see, e.g. BT87, p. 205; Tonry, 1983; Merrifield and 
Kent, 1990; Kent, 1990). Indeed using the tomography the calculations are done 
using only cartesian coordinates. In our particular case, due to the linearity of the 
with respect to the angular coefficients and velocity monomials [see (12)], it is 

evident that: 

[v^]n = [v^]n. (19) 

So we can equally use or in any discussion involving the angular mean 
operator. In order to evaluate (19) for iV > 0, we need the following 

Lemma For any integers N > K > L > then: 




N-K K-LLi 

"i "2Jn 



T[{N + 2)/2][l + (-1)^-^][1 + (-1)^-^][1 + (-1)^] 
8(iV + l)T[{N -K + 2)/2] T[{K -L + 2)/2] T[{L + 2)/2] 

(20) 



Proof. The angular mean is calculated using formula 8.380.2 of Gradstheyn 
and Ryzhik (1965). Then, using (^) = T{a + l)/[T{b + l)T{a +6-1)] and the 
duplication formula for T, after some algebra the result is proved. 

This gives the following: 

Theorem For any integer iV = 2n + 1 > 0, 

[v^'h = 0; (21) 

for any integer N = 2n > 0, 



1 ^^fn\(k 

[-""In 



"'--.^^^EE (.rV-'4)^ (22) 

fc=o;=o \ / \ / 

Proof. From (12) it follows that in each addendum of the sum of the 
exponents is equal to N, and from (20) only the terms with three even exponents 
have angular mean different from zero. For iV = 2n + 1 at least one among N — K, 
K — L and L is odd. For N = 2n the numbers N — K, K — L, and L, are necessarily 
either all even, or one even and two odd: setting N = 2n, K = 2k, L = 21, also 
(22) is proved. 

For example, for N = 2, [v'^]q = {vf + t^f + i'|)/3, and from (11) 



[^2]^ = [^2]^ = - ' - ' - + li^^:^^^. (23) 
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In particular the identity [(T^Jq = [af^ + (tIj + (t^^)/3 is proved. The analogue of 
(15), i.e. the equation for the angular mean value [t;^]n is immediately obtained: 



dt dxi V dvi ) dxi 

5 Tomography of the general PVT 

We can now obtain the general expression of the angular mean value of the PVT 
of any order for stellar systems of unspecified shape from (16). The angular mean 
value of the projected kinetic energy is obtained from (14) using (21) and (22). 
The mean value of the generalized gravitational energy in (25) is apparently more 
complicate. On the contrary it is immediately obtained from (17) observing that 



1 d^[v^+^] 



and using (21) and (22). For example, the tomography of the classical PVT is 
given setting iV = 1 in (16): 

W 

2[K2,iosh = -[Wo,ios]n = -y. (26) 

Considering now a spherically symmetric system the integral of the l.o.s. square 
velocity weighted on I divided by C is equal to 2/3 of the total kinetic energy of the 
system; moreover, if the system is characterized by pure velocity dispersion, (1) is 
naturally recovered. 
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